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We study a two-state statistical process with a non-Poisson distribution of sojourn times. In 
accordance witii earlier work, we find that this process is characterized by aging and we study three 
different ways to define the correlation function of arbitrary age of the corresponding dichotomous 
fluctuation. These three methods yields exact expressions, thus coinciding with the recent result by 
Godreche and Luck [J. Stat. Phys. f04, 489 (200i)]. The first method rests on the lines established 
in an earlier paper [AUegrini et al., Phys. Rev. E 68, 056f23]. These authors built up an infinitely 
aged Generalized Master Equation (GME), compatible with the ordinary Onsager Principle. Herein 
we generate first a GME of arbitrary age, then we derive from it an expression for the corresponding 
correlation function, based on the Onsager principle extended to a condition of generic age. To derive 
the second exact expression of correlation function of arbitrary age, we use a Liouville-like approach, 
based on a model mimicking the environment responsible for the fluctuations of the dichotomous 
variable under study, and producing the same non-Poisson distribution of sojourn times as that 
considered in the earlier treatment. This exact expression refers to time t, while the earlier exact 
expression refers to the variable u, which is the Laplace conjugate of t. Finally, the third exact 
expression, again with respect to time, is derived from the direct observation of the dichotomous 
sequence generated by the two-state system fluctuations. The resulting expression is implicit, and 
we adopt a convenient approximation to obtain a simple analytical formula. This approximation 
rests on the renewal nature of the process, which resets to zero the memory of the system after 
the occurrence of any event, thereby suggesting that a Markovian approximation can be made. 
Actually, non-Poisson statistics yields inflnite memory at the probability level, thereby breaking 
any form of Markovian approximation, including the one adopted herein. For this reason, we check 
the accuracy of the analytical formula by comparing it with the numerical treatment of the second 
of the three exact expressions. We flnd that, although not exact, the simple analytical expression 
for the correlation function of arbitrary age is very accurate. 

PACS numbers: 02.50.Ey, 05.40.Fb, 05.20.-y, 05.60.Cd 



I. INTRODUCTION 



The phenomenon of aging has been known for a long 
time as a property of spin glasses and polymers 1]. Part 
of the reason for the more recent interest in this phe- 
nomenon has to do with the predicted breakdown of cer- 
tain fundamental assumptions made in equilibrium sta- 
tistical mechanics when applied to strongly disordered 
systems. For example, the Onsager Principle 0, being 
the relaxation of a perturbed system back to its equi- 
librium state as described by an unperturbed autocor- 
relation function, is violated in anomalous diffusion and 
anomalous relaxation. More recent papers on this phe- 
nomenon are devoted to studying aging in diffusion pro- 
cesses occuring in d-dimensional lattices in low di- 
mensional environments ^ and in the quantum dynam- 
ics of dissipative free particles . Most recently 0, |7| 
there has been some interest in the manifestation of ag- 
ing in processes described by means of the Continuous 
Time Random Walk (CTRW) formalism Q. 

A natural way of formally expressing aging is through 



the correlation function of the stochastic variable of in- 
terest, ^, which is expressed as an ensemble average in- 
dicated by the brackets 



t2). 



(1) 



The averaging brackets carry a subscript stressing that 
this is a ^a-old property, rather than the traditional aged 
property, only depending on \ti ~ t2\. The time differ- 
ence indicates the stationarity of the underlying process, 
but the ta-superscript denotes a dependence on the age 
of the process as well. In the case of a dichotomous 
variable obeying renewal theory, the exact expression for 
the Laplace transform of this age-dependent correlation 
function was found by Godreche and Luck "o*]. Their ex- 
act result was recently recovered by Margolin and Barkai 
|lfl| | as a special case of a more general expression, since 
the latter authors do not require the condition that the 
two states of the dichotomous variable ^ have the same 
waiting-time distribution. However, herein we make the 
same assumption as did Godreche and Luck 0, and de- 
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rive their result alongthe lines of the recent work of Ref. 
0. AUe erini et al. |3 noticed that in the non-Poisson 
case the well known Generalized Master Equation (GME) 
of Kenkre, MontroU and Shlesinger fllj becomes incom- 
patible with the Onsager principle (21 ^^id found a way 
to make the GME compatible with the aged condition. 
However, they 01 l^ft unsolved the problem of deriving 
a GME of arbitrary age, which is equivalent to making 
the Onsager principle compatible with an incompletely 
aged system. Here we solve this problem and establish 
that this solution leads to the exact expression for the 
correlation function of arbitrary age found by Godreche 
and Luck 

The GME derived herein does not have the same ori- 
gin as those widely discussed in the literature 0. The 
Zwanzig GME, in fact, one of the most popular mas- 
ter equations, is derived from a first principles proce- 
dure, starting from the statistical Liouville equation of 
the whole universe. The Zwanzig GME is a projection of 
this universal Liouville equation onto the Hilbert space of 
the system of interest. The approach used to derive the 
GME herein is quite different from that of Zwanzig and 
is based on the experimental observation of non-Poisson 
dichotomous signals. Examples of such signals include 
those produced by an ionic channel _13j and by blinking 
quantum dots We build up a GME that is com- 

patible with the experimentally determined non-Poisson 
nature of these processes, assuming the applicability of 
renewal theory. We leave open the question as to the 
source of randomness, but note in passing that the fluc- 
tuations in the system variable are generated by the envi- 
ronment. A simplified dynamical model has reproduced 
the essential statistical properties of the system of inter- 
est, in the case of infinitely aged GME In the case of 
arbitrary age, not only a first principle derivation is miss- 
ing, a condition shared by the infinitely aged case: there 
is, to date, not even a heuristic derivation derivation of 
the corresponding GME. Herein we provide a heuristic 
derivation of a GME with arbitrary age, one based on 
an exact expression for the corresponding distribution of 
sojourn time, and we use this result to derive an exact 
expression for the correlation function. 

In conclusion, we point out that we rework the prob- 
lem from a number of perspectives. We do this because 
each approach provides separate and distinct insights into 
the phenomenon described by a GME of arbitrary age. 
The first perspective adopted in this paper focuses on the 
derivation of a GME, based on the experimental obser- 
vation of the time evolution of a trajectory, characterized 
by rare jumps from one state to another. The second per- 
spective uses a Liouville-like approach to the time evolu- 
tion of the variable of interest and its environment. No 
attempt is made to derive the former dynamic picture 
(trajectory) from the latter (probability density). The 
probability density perspective is characterized by infi- 
nite memory, yet, the statistical process under study is 
generated by random critical events, whose occurrence 
erases the memory of earlier events. Therefore, we think 



it is prudent to examine, yet again, the same process from 
a third perspective, one based on the direct observation 
of the sequence of rare random events. 

We find that the use of three different perspectives is 
fruitful in shedding light on the recent observation made 
by Sokolov, Blumen and Klafter Their results, in our 
opinion, imply the breakdown of certain well-established 
notions of equilibrium statistical mechanics, such as lin- 
ear response theory as a prescription for predicting the 
effect of externally perturbing a system out of equilib- 
rium. 

The outline of the paper is as follows. In Section II we 
employ the reduced density perspective, namely, we build 
up the GME of arbitrary age. In Section III we study the 
time evolution of the total distribution density. In Sec- 
tion IV we use a third method, based on the probability 
of occurrence of the rare random events. Section V is 
devoted to concluding remarks. 



II. THE REDUCED DENSITY PERSPECTIVE 

This section is devoted to the application of the re- 
duced density perspective to the construction of the 
GME. In the literature of statistical physics this per- 
spective implies a contraction on the total density onto 
a prescribed subspace. Consequently, constructing the 
equation of motion of the total density matrix should 
be the first step. The reason why we reverse the per- 
spective, thereby confining the total density treatment 
to the next section is given by the fact that non-Poisson 
statistics has the effect of making the ordinary approach 
extremely difficult, if not impossible, as discussed in the 
companion papers [11 113. The derivation of the GME 
is made possible by expressing the higher-order correla- 
tion functions in terms of the second-order correlation 
function, via simple expressions that are violated by the 
Poisson statistics. This is the reason why we derive the 
GME from the Continuous Time Random Walk (CTRW) 
of MontroU and Weiss rather than from the equa- 

tion of motion for the total distribution density. 



A. From the generalized master equation to the 
age-dependent correlation function 

The reduced density matrix perspective is based on the 
adoption of the GME formally defined by 

j^p{t)^- j\{t-t')Kp{t')dt', (2) 

where p (t) is the m-dimensional population vector of 
m sites, so that its ith component, Pi{t), represents the 
probability of finding the walker at time t in the ith site. 
K is a transition matrix between the sites and $(t) is the 
memory kernel. We study the case where the fluctuating 
variable ^ leading to this GME is dichotomous. This 
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means that the GME has only two states and the matrix 
K reads 



K = 



1 -1 
-1 1 



(3) 



The waiting-time distribution in either of the two states is 
denoted by ■ip{t). We estabhsh the connection of the GME 
with the CTRW 0i 111 by relating the Laplace transform 
of the memory kernel given by 't{u), to the Laplace 
transform of the waiting-time distribution ')jj(t), given by 
as follows 



$(w) 



uip{u) 
l-ip{u)' 



(4) 



We consider two distinct processes, that of prepara- 
tion and that of observation. The preparation process 
establishes the initial conditions of the set of trajectories 
under study. A trajectory is a sequence of symbols -I- 
or — , specifying whether the system is in the state |-|-) 
or |— ). We call the time interval with the system ei- 
ther entirely in the state |-|-) or entirely in the state |— ), 
the laminar region. With Zumofen and Klafter we 
assume that the preparation process, beginning at time 
t = —ta, insures that all these trajectories begin with the 
system at the onset of a laminar region, either + or — . 

The observation process begins at time t = > —ta- 
The distribution of first sojourn times is denoted by 
iptai't) ^^'^ is, in general, different from V'(i)- In fact, the 
laminar region corresponding to the first sojourn time 
might have begun earlier than t = 0. The exact expres- 
sion for this time distribution is 

i^tAt) = + ta) + Yl / dyij{y + t)iJn{ta - v), (5) 

where ipniT) denotes the probability that n jumps occur 
during the time interval of length r, the last of which 
occurs at time t — t. As is well known, for a renewal 
process the waiting times for successively more jumps is 
given by the convolution 



i>n~l{t')iji{t - t')dt' , 



(6) 



with V'(i) = V'l (0- The io-old distribution of first sojourn 
time was discussed earlier in some detail 0. However, a 
careful analysis of Eq. (5) can help the reader to realize 
the rationale behind this distribution of sojourn times. 
The first term on the right hand side of Eq. |(SJ corre- 
sponds to the case when the first laminar region is ex- 
tended in time more than the time interval ta between 
preparation and observation. The second term takes into 
account the cases when the last laminar region began af- 
ter the preparation time, after a sequence of an arbitrar- 
ily large number of earlier laminar regions, the first of 
which, of course, begins ai t = —ta- 



in the Poisson case, because of its unique functional 
form, there is no dependence of il^t^ on ta, and conse- 
quently no aging. In the non-Poisson case, on the con- 
trary, the two waiting-time distributions, 
are identical only iita — Q- In this case both CTRW and 
GME correspond to switching-on the observation process 
at the same time as the preparation process, and the con- 
nection between the two pictures is given by Eq. |0J. 

AUegrini et al. proved that the GME is compatible 
with an infinitely aged CTRW, provided that the memory 
kernel is made compatible with an infinitely aged 
condition, characterized by a distribution of first sojourn 
times, which is infinitely aged. In this case ia — *■ oo so 
that 



*oo(u) 



l+i^iu) - 21/^00 (w) 



where 



1 r°° 

^oo{t) = / dt'^j{t' 



(7) 



(8) 



It is straigthforward to extend the calculations of Ref. Q 
to the case of an arbitrarily ^a-old system, so that Eq. Q 
is replaced with 



(9) 



where tptai^) is the Laplace transform of ipt^{t). 

Using the GME with the ia-old memory kernel, we 
define the age-dependent correlation function <I>^*''^(t) 
through its Laplace transform, as follows 



1 



u + 2$t„(u) 



This prescription corresponds to setting 



Plit)-P2it) 
Pl(0)-P2(0)' 



(10) 



(11) 



as one can easily check by Laplace transforming both 
sides and using the GME with the ta-o\d memory kernel 
(see section II B). In other words, the correlation func- 
tion of arbitrary age mirrors the extension of the Onsager 
principle, which is usually limited to infinitely aged sys- 
tems 1^, to physical conditions of arbitrary age. It has 
to be pointed out that at this stage there is no guarantee 
that the Onsager principle of arbitrary age holds true. 
However, in Section IIIBI we show that Eq. (|10|l yields 
the exact result of Godreche and Luck ^] for the corre- 
sponding correlation function. Furthermore, in Section 
III CI we establish the same result using arguments based 
on trajectories rather than densities, thereby affording 
an independent construction of the exact expression for 
the ta-o\d correlation function. All this can be thought 
of as a compelling demonstration of the correctness of of 
Eq. extending the Onsager principle to conditions 

of arbitrary age. 
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B. Derivation of the exact expression proposed by 
Godreche and Luck: the probability perspective 

To establish that the proposed approach yields the ex- 
act expression of Godreche and Luck let us express 
the ta-o\d correlation function through the probability 
vector p = (pi,p2) for the dichotomous variable f = ±1 
to have either positive (state 1) or negative (state 2) val- 
ues. We have that: 



mm))t 



(12) 



Pi(0)pi(i|l,< = Q)+p2{Q)p2{t\2,t^ 0) 
-Pi(0)p2(i|l, t = 0)- p2(0)pi(i|2, t = 0), 



wherepj(i|fc, t = 0) is the conditional probability that the 
variable ^ is in the state j at time i, given that at time 
t = it was in the state k. This means that pj(t\k, t — 
0) is obtained letting those trajectories evolve that at 
time t — had ^ in the state k. For a straightforward 
evaluation of pj{t\k,t = 0), we use the GME formalism, 
adapted to the ia-old system, and we take into account 
the initial condition pi(0|fc,t = 0) = 5i^k- According to 
the GME the components of the conditional probability 
vector are determined by: 

d /■* ^ 

-p^{t\k,t = 0) = - / dt'<^t^{t~t')Y,Kj,p^{t'\k,t = 0), 

•^0 i=l 

(13) 

with {u) given by Eq. Q and the elements of K given 
by Eq. ®. 

By Laplace transforming (|13|l and doing some algebra, 
we obtain 



p,(7/|fc,0) =^(wl-f <lt„(u)K)T/p,(0|A;,0). (14) 
1=1 

Defining the matrix 

( "+'i'to(") ^ta{u) 
^tg (m) "+^*a (") 

u[jj+2<i>t„ («)] u[«+2<it„(u)] 

and using the initial condition pi{Q\k,t = 0) = 5i^k, we 
obtain for the Laplace transform of the conditional prob- 
ability vector the following expression 



Pj{u\k,Q) 

(15) 

Using Eq. (|15|l we can Laplace trasform (|12|) to obtain 



M + $t„('u) 




k + ^t„(u) [5j,k+i + Sj^k-i) 


U 







e (")=pi(o) 



u (u + 2^tAu) 



'Pi(O)- 



u(u + 2$tju) 



-P2(0) 

-P2(0) 



u(ii + 2<I>f„(u) 



u{u + 2^tAu) 
(16) 



We note that the probability is normalized, pi{Q) 
P2{0) = 1. Thus, it follows that 



u + $t„(u) 



+P2{ 



u(u + 2$t^(u)) u(u + 2$t^(u) 
1 



M-l- 2$i^(w) ' 



bi(0) 



(17) 



confirming the correctness of the definition introduced in 
Eq. Substituting ^ into we obtain 



$[*°^(ii) 



1 



1 



-^tg (") 



2iJtAu) 

1 + tpiu) 



(18) 

which coincides with the results of Godreche and Luck . 
Furthermore, the ratio of the differences in probability is 
determined in Laplace space by 



Piju) -P2{u) 

Pl(0)-P2(0) 



{Ju - ■hi)Pi{^) 



1 

u + 2$f„(u) 



P2(0) 



(19) 



= l>J*"'(u). 



As pointed out earlier, Ea. (|19|l means that one can ex- 
tend the Onsager Principle from the infinitely aged sys- 
tems, for which Onsager originally defined it, to systems 
of any age. In the latter case the relaxation is propor- 
tional to the ta-old correlation function, not to the in- 
finitely old, or equilibrium, correlation function. In sum- 
mary, we discovered an Onsager principle of arbitrary 
age, at least in the special case of the dichotomous vari- 
ables considered in this paper. 



C. Derivation of the exact expression proposed by 
Godreche and Luck: the trajectory perspective 

It is possible to again derive the exact result of Eq. 
(|18|l from a different perspective, which will allow us, 
in Section llVI to propose an analytic expression for the 
ta-old correlation function as a function of time. This 
expression, as we shall see, is not exact, but it is shown 
numerically to be a very good approximation to the exact 
result. 

The usual method of connecting the correlation func- 
tion $j to the waiting-time distribution, within a tra- 
jectory perspective, is to introduce a theoretical waiting- 
time distribution, ip*(t), which cannot be observed di- 
rectly. In fact, the experimental waiting-time distribu- 
tion, namely, the distribution of times with alternate 
signs, denoted by us as ?A(i), is obtained from the theo- 
retical waiting-time distribution, ^l^*{t)^ by adopting the 
following procedure. We divide the time axis into bins, 
whose size is determined by the waiting-time distribution 
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tp*{t). Then, these bins are assigned either the value 1 or 
the value —1, by tossing a coin to make the decision. It 
is evident that the intervals along the time axis with the 
same sign, are larger than the time bins determined by 
'4'*it), since two or more consecutive coin tossings might 
have produced the same sign. It is shown that the 
Laplace transform of tJj is connected to the Laplace trans- 
form of ■(/'* via the relation 



(20) 



Let us use the term event to denote the coin tossing 
introduced above. The expression H20|l is the result of 
summing over all possibilities of not changing sign with 
a coin toss, which turns out to be a geometrical series in 
the Laplace representation. The correlation function <I>^ 
and the theoretical waiting-time distribution function ■(/;* 
are connected through the relation 



{r) 



where the average waiting time is given by 



(t) — I T1p*{T)dT. 

Jq 



(21) 



(22) 



Eq. H21|) determines that the correlation function ^^{t) 
is equal to the probability of finding a window of length 
t without internal events. 

The same result can be immediately recovered using 



where the symbol * denotes time convolution. In fact, 
after a first interval of time followed by a coin toss with no 
change of sign, determined by V't„/2, the next intervals of 
time with no change of sign according to the coin tossing 
prescription, are determined by ip* /2. The sum of the 
convolutions takes into account all the possible sequences 
of intervals of time with no change of sign before a change 
of sign of the variable ^ eventually occurs, and gives as 
final result the distribution for a first observed sojourn 
time t of the variable ^ in one of its two states, that is 

Thus by summing the geometric series in the Laplace 
variables, from H26|l we obtain 



2 - ^*{u) 
Using Eq. if^ we write Eq. (|77ll as 



i + '4'iu)' 



(27) 



(28) 



By Laplace transforming Eq. H25|) and using Eq. H28|l . 
we obtain 



d[*"'(u) 



l + ^JJ{u) 



(29) 



namely, we again recover the exact result of Godreche 
and Luck given by Eq. (|18|l . This establishes the equiv- 
alence of the trajectory and GME prescriptions for this 
process. 



1 pOC nOO 

^ TT dt' dt"r{t")- 
{t) Jt Jt' 



(23) 



We note that, see Ref. 0, 

rAt)^yJ^ dt'rit') (24) 

is the infinitely aged waiting-time distribution. Actually, 
Eq. H23|l is the infinitely aged correlation function, a 
special case of the more general prescription 

/oo />t 
dt'rtSt') = i- dt'rtSt'). 

(25) 

It is straightforward to show that the ta-old experi- 
mental waiting-time distribution and the ta-o\A theoret- 
ical waiting-time distribution are connected through the 
following sum of convolutions: 



(26) 



D. Generalized Master Equation of arbitrary age 

We are now in a position to make a preliminary balance 
of the results obtained so far. The first is that we have 
generalized the result of an earlier paper 0, in that we 
have derived the GME of arbitrary age 



d , ^ 



'^tAt~t')¥.^{t')dt', 



(30) 



whose memory kernel [t) is defined though its Laplace 
transform by means of Eq. @ . We have also shown that 
the Onsager principle, valid for infinitely aged systems, 
can be extended to conditions of any age, and that this 
extension allows us to derive an exact expression for the 
ta-old correlation function. However, the analytic results 
obtained so far are in the Laplace domain. It is desireable 
to achieve them in the time domain, as well. We now turn 
our attention to the latter. 



III. THE LIOUVILLE-LIKE APPROACH 

In this section we derive another expression for the 
ta-old correlation function adopting a perspective where 
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aging is determined by the out of equilibrium bath for the 
variable of interest, ^. There might exist conditions, as 
we shall see, where equilibrium is not even allowed. The 
expression for the ia-old correlation function afforded by 
this perspective is exact, and is thus equivalent to the 
Godreche and Luck expression of Eq. H18|l . However, the 
exact expression is implicit, rather than explicit, and is 
therefore more convenient for the numerical calculations 
done subsequently. 

Here we adopt the perspective of earlier work |0, |^ 
to account for the aging effects characterizing the fluctu- 
ations of the dichotomous variable ^. These fluctuations 
occur while the environment of the variable, ^, slowly 
drifts. This drifting process is extended over time, and 
could lead to circumstances where it is not possible to 
attain equilibrium asymptotically. In keeping with the 
jargon of statistical mechanics the environmental or "ir- 
relevant" variable is called y and in the model moves in 
the interval / = [0,2]. In the semi-interval [0, 1], we use 
the equation of motion for the probability density 



(31) 



This is the motion determined by a potential, with the 
minimum at y = 1, in the over-damped case. In the 
interval [1,2], the over-damped potential is the mirror 
image of the potential acting on the left interval. Conse- 
quently, if the initial condition is located in the internal 
part of the interval [0,1], the particle moves, from the 
left to the right, with a deterministic motion, towards 
y = \. If the particle is initially located in the interior 
of the interval [1,2], it moves deterministically from the 
right to the left. When the particle reaches the poten- 
tial minimum, it is injected back, with equal probability, 
into any of the points of the interval / excluding y = 0, 
y = \ and y ~ 2. The time spent by y within / cor- 
responds to sojourning in one of the two states of the 
variable ^, either |-|-) or |— ). The instant of back in- 
jection corresponds to the choice of the new state and, 
with equal probability, this is either the same state or the 
other state. The variable y represents the environment 
of the variable ^, and its initial distribution is given by 
the state of the bath. The corresponding waiting-time 
distribution between two consecutive back injections is 



rpfj,- 



where the index ^ is related to z of Eq. H31() by 



z-1' 



(32) 



(33) 



and the parameter T, characterizing the waiting-time dis- 
tribution of Eq. is 



T ; 



A^-1 
A ' 



(34) 



in accordance with the normalization constraint. The au- 
thors of Refs. f^jlsl have shown that the essential prop- 
erties of the dichotomous non-Poisson fluctuation can be 



accounted for by limiting ourselves to this simplified pic- 
ture, involving only the semi-interval [0, 1]. In this simple 
picture the aging process is described by 

^^ = -A^2/M2/,i) + Ap(l,t), (35) 

which takes into account the back injection into the semi- 
interval, occurring with uniform probability, when the 
particle reaches the point y — I. For simplicity, but with 
no loss of generality, we fix A = 1 . 

Using the results of Ref. jlTj the ta-old correlation 
function is evaluated as follows. The bath is prepared 
at time t = ~ta- This means that at time t = —ta the 
distribution of y within the interval / is flat. In the case 
z < 2 this distribution tends to the equilibrium distri- 



bution, Pe 



If z > 2 this distribution diverges. 



thereby implying that the distribution approaches the 
Dirac delta function located at y = 0. This is the non- 
stationary condition, the condition where equilibrium is 
not allowed, and only a condition of eternal drift is ad- 
mitted. Suppose this distribution evolve for a time ta, 
without our observing it. This means that we begin the 
observation process when the system has a new distri- 
bution, different from the initial flat distribution, and 
determined by its time evolution from t = —ta to t — Q, 
described by Eq. For i > Eq. (jSHl) is replaced 

with Eq. ijni) 



dp{y, t) d 



dt 



dy 



(36) 



namely the back injection process is stopped, thereby 
implying that the population decreases. 

The theory of Ref. 0] relates the probability solu- 
tion to the Liouville-like equation to the ia-old correla- 
tion function as follows 



$(*"'(t)= / dyp{y,t). 



(37) 



Note that the initial condition p{y, 0) is obtained from 
Eq. H35|l moving from the flat distribution at t = —ta- 
The time evolution, corresponding to the observation 
process, is determined by Eq. (|36() . This simple pic- 
ture is a fair representation of the description made in 
terms of trajectories in Section |nj The fact that the 
norm of p{y,t), when p{y,t) is described by Eq. (|36|l . is 
not conserved, reflects the occurrence of jumps from one 
state to the other, making the population of a given state 
decrease. Let us remark that the correlation function is 
determined by the antisymmetric part of the whole dis- 
tribution |l7j ]. This formal condition corresponds to ob- 
serving the time evolution oip{y, t), with y ranging from 
y — Oto y — 1, under the action of Eq.||31J, with no back 
injection process. The process of back injection is es- 
sential to determine the correlation function of arbitrary 
age. However, once the parameter ta is fixed, and with it 
the age of the correlation function, the evaluation of the 
correlation function is done by imagining the bath frozen 
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in the distribution corresponding to this fixed age. This 
is the reason we use Eq. (|31f) to determine the correlation 
function rather than Eq. 1)35(1 . 

The process of back injection is essential for the slow 
environmental drift, but the correlation time of age ta 
is determined by the age-fixed condition, the bath being 
id-old, and keeping this age forever. 

Our goal is to calculate the correlation function of arbi- 
trary age. According to Eq. ((37(1 . this goal requires that 
we evaluate p{y,t), first. To carry out this calculation 
we follow the procedure illustrated in detail in Ref . . 
The time evolution from t = —ta to t < is described by 
Eq. ((351) . yielding for the solution 



1 



[l + (^-l)(t + ta)2/^-l]^/(^-l) 

* p(.l,r)dT 

[l + (z-l)(t-r)j/^-i]^/(^-i)' 



(38) 



with p{l, t) determined, in turn, by Eq. ((35() . 

As already stated, for t > 0, p{y, t) is governed by Eq. 
((36(1 . so that its solution is, given the initial value. 



[l + {z- l)y^-it]^/(^-i) "'-[i+(.-i)y--it]i/(--i) 

(39) 

Now, we can calculate the final expression for the age- 
dependent correlation function given by Eq. (|^ . After 
some simple algebra we obtain the following expression: 



1 



f" p{l,r)dT 
J-U [l + (z-l)(t-r)]i/(-i) 

^ P{l,T)dT 



(40) 



[l + (^-l)(<-r)]i/(-i) 

^*{t + ta) + [l-^*{ta)] 

" p(l,r)dT 

[l + (z-l)(i-r)]i/(-i) 

p(l,T)dT 



[l-(z-l)r]i/(-i)- 
Note that in Ea. l(4l)() we are using the definition 



• if the age vanishes — > 0, then the ta-o\<i correla- 
tion function reduces to the probability of no event 
**(t). 



occuring <bf''\t) 



• if the age increases without limit ta ^ oo, then 
the ta-old correlation function reduces to a known 



result $l*°-'(t) ^^(t), where 



[l-h(z-l)t](2-^)/(^-l) 



(42) 



is the usual correlation function calculated at equilib- 
rium. So, we must observe a crossover between two dis- 
tinct regimes. Moreover, we anticipate that we recover 
these results from a different perspective in Sec. IIVI 

The perspective adopted in this section makes it eas- 
ier to understand, on physical grounds, why the Poisson 
condition produces aging. The Poisson condition corre- 
sponds to z = 1 in the nonlinear dynamical map, and 
in this case the distribution is always fiat. Aging is the 
process of slow regression of the bath variable to equilib- 
rium, a process that can, in principle, last forever. With 
this representation it is possible to study the correlation 
function of arbitrary age also in the case p, < 2. AUe- 
grini et al. (ITl] give a detailed discussion of the effect of 
external perturbations. 

In the case z < 2 equilibrium is possible. An external 
perturbation can be used to create a condition whereby 
the system is out of equilbrium. The observation of the 
population difference between the two states, ^ = 1 and 
^ = — 1, is equivalent to determining the correlation func- 
tion. Note that the initial condition that we use to study 
the time evolution of p(y, t) during the observation pro- 
cess, Eq. (P5|) . is the antisymmetric distribution of the 
interval [0,2] folded into [0,1] (see Ref. Tt^). It is ev- 
ident that in practice the equilibrium correlation func- 
tion is never observed. In fact, the external perturbation 
should create a non-symmetric initial condition in the 
whole interval [0,2]. The anti-symmetric portion of this 
distribution, folded into the interval [0,1], should create 
the same out of equilibrium distribution as that produced 
by folding in the same way the anti-symmetric portion of 
the equilibrium distribution. As pointed out in Ref. |l7j . 
this is impossible to do in practice with a a perturba- 
tion lasting for a finite time, and, thus, with any realistic 
experimental procedure. 



5-*(t) = j ^*{T)dT. 



(41) 



This notation for the survival probability is also used in 
Section FlV Al We notice that, t > and p{l,t) is calcu- 
lated from Eq. (j^HJl with y ^ 1 and -ta < t < 0. This 
is an exact expression that turns out to be convenient to 
numerically check the prescriptions of Section IIVI 

It is worth using the perspective of this section to shed 
further light into the problem under discussion in this 
paper. Let us notice that 



IV. AN ANALYTICAL EXPRESSION OF THE 
CORRELATION FUNCTION OF ARBITRARY 
AGE 

The exact expressions, Eq. ((T^ and Eq. al- 
low us to determine the dependence of the age-dependent 
correlation function on time via anti-Laplace transforms. 
This inversion is done numerically since we are unable to 
carry out direct analytic inversion of the Laplace expres- 
sion. There is therefore the need to find an analytical 
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expression for even if the price we pay is an ap- 

proximation and the loss of exactness. This section is 
devoted to the derivation of an accurate, albeit approx- 
imate, analytical expression for the correlation function 
of arbitrary age. 



An exact implicit expression for the correlation 
function of arbitrary age 



We write namely the probability of not finding 

events inside the interval (0, t), given the occurance of an 
event at time —ta, in the following implicit form 



dtir{ta-tl)<^f\t). (43) 



This is the probability of finding an event at a time larger 
than t from the preceding one. The right-hand side of 
Eq. (|43|l is the sum of two terms, the first corresponds 
to no event occurring in the whole interval {—ta,t) and 
the second term takes into account the possibility of at 
least one event occurring between —ta and 0. The instant 
—ti signals the first of these events (or the event, if there 
is only one). The coin-tossing procedure has the double 
effect of rejuvenating the process (the correlation func- 
tion now has age U) and of factoring the two probability 
functions inside the integral. 

It is interesting to notice that Eq. (|43|l is exact, which 
can be verified by differentiating both sides of the equa- 
tion, and making use of H25|l and H41(l . This procedure 
yields 



rtAt)^r{t+ta)+ dhrita-hWtAt), (44) 

Jo 

namely an implicit (but analytic) expression for . A 
straighforward iteration of (|44|l leads to 

rtM = r{t+ta)+ r dhr{ta-ti)r{t + h) 

Jo 

+ f ° dhrita ~ti) r dt2r{tl - t2W{t + t2) 

Jo Jo 

+ f ° dhrita -h) r dt2r{h - h) 

Jo Jo 

X r dhr(t2-h)r(t + h) + -- - , 

(45) 

i.e., a sum over all possibilities of finding events at both 
times t and —ta, with no events between and t. This is 
exactly the definition of lAt* [t]- In detail, the first term 
in the right-hand side of corresponds to the case in 
which no event occurs in the interval / = {—ta, 0], while 
all further terms refer, respectively, to one event at ti 
in /, to two events, at ti and t2, in /, and so on, with 
ti < t2 < t^ < ■ ■ ■ . Note that ipl^ [t) is obtained from 



t]-'*{t) in the same way iptait) is obtained from ip(t): in 
fact if in Eq. (|3J) we make the substitution ip{t) i'*{t) 
we obtain exactly Eq. H45|) . i.e. the definition oiip^^{t). 



B. Equivalence with the exact expression of 
Godreche and Luck 



As we said, Eq. is exact. Consequently, this equa- 
tion is expected to be equivalent to the exact expression 
proposed by Godreche and Luck in Ref . [7] . In this sub- 
section, as a double check, we prove that Eq. (|43|l is in 
fact equivalent to the expression of Godreche and Luck. 
To do that, let us differentiate both sides of this equation 
with respect to t to obtain 



|$(*"'(i) = -V*(i + U 



dtirita-tl)j^^'^^'\t). 



(46) 

A more tractable expression is obtained by taking the 
Laplace transform of this new expression with respect to 
t„ to obtain 



^*(s) 



dye'^yriy) 



(47) 



+ r{s)j^^'\t), 



from which it follows that 

e^'\^>*{s)-J^dye~^yr{y) 



1 - V*(s) 



(48) 



The right-hand side of this equation appears formidable, 
but one can easily see by Laplace transforming Eq. 144|l 
with respect to ta, that it is nothing more than the 
Laplace transform with respect to ta of — V't (^)i namely 
—^p*(t). Consequently, by inverse Laplace transforming 
Eq.(||HI) we obtain 



-^['^\t) = -rtAt), 



which, by integration, yields 



$f)(t) (0)- / dt^^Ut,) 



(49) 



(50) 



The initial value of the ta-old correlation function can be 
determined by explicitly Laplace transforming Eq. H43|l . 
calculated at t = 0, with respect to ta, thereby yielding 



which simplifies to 

1 - V^*(s)j $['^(0) = (s). 



(51) 
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On the other hand, we know that ^* (s) 
1 — 7/>*(s) /s, so the initial Laplace variable is 



1/s. 



It is obvious that the inverse Laplace transform of this 
last expression yields the initial condition for the corre- 
lation function $^*''''(0) = 1, as it must be, based on the 
definition of Eq. when ti — t2 and ^ = ±1, a di- 
chotomous variable. Using this initial condition in (|50|l . 
it follows that 



Jo 



(52) 



which, as we know from Section II, coincides with the 
exact expression for the ^a-old correlation function of Ref. 

i 



This expression is not exact, but it is analytic. Moreover, 
the form of Eq. suggests an iterative approach, and 
we can therefore refine our result, by replacing <I>^' (t) with 
A^*'^\t), after an arbitrary number of iterations. In Fig. 
1 we check the accuracy of the first-order approximation 
to the exact expression given by (|54|l . by comparing it 
to the exact prediction of Eq. H40() . The curves in Fig. 
1 represent a numerical treatment of 1)40(1 with different 
values of ta- In this example = 2.5 and T = 1.5 (as 
said, A = 1). The highest curve in the figure represents 
the stationary correlation function, namely ta = oo. The 
correlation function $^*"''(t) has a faster decay, with de- 
creasing values of ta ■ Finally, the lowest curve represents 
the correlation function with zero age, namely ^'^ (t). 
The symbols in Fig. 1 overlaying the continuous curves 
represent the calculations using Eq. H54II with various 
ages. We see perfect agreement for ta < (t) and for 
ta ^ (t)- The agreement remains good for intermediate 
values, with discrepancies comparable to the numerical 
round-off errors. 



Approximation through iterative expansion, 
and truncation 



Although we established that Eq. 1(43(1 is exact, due 
to its implicit nature, it cannot be used to obtain an 
analytic function without carrying out a reasonable ap- 
proximation. The approximation we select rests on re- 



placing <i>^*'^(t) in the right-hand side of Eq. 

the function yl'^*"-' (t), which denotes a correlation func- 
tion of uncertain age, which is, however, younger than 
the ta-old correlation function defined by Eq. 1(43(1 . This 
correlation function reads 



with 



(53) 



The easiest way to derive ((53(1 is to adopt the language of 
conditional probabilities. The numerator of this expres- 
sion is the probability of not finding an event between 
and t, minus the probability of not finding an event be- 
tween —ta and t. If we call A the condition of no event 
in the interval (0, t), and B the condition of at least one 
event in the interval {—ta,0], then the numerator of the 
right-hand side of ((53(1 can be identified with the joint 
probability P{A,B). On the other hand, the denomina- 
tor of 1(53(1 is simply P{B), the probability of the event 
B occuring, and therefore the function A(*")(t) can be 
identified, as it should, with the conditional probabilily 
P{A\B), namely with the probability of finding no event 
between and t, {i.e. it is a correlation function), given 
an event between —ta and {i.e. it has an age younger 
than ta)- 

Finally, by plugging Eq. ((53(1 into Eq. 1(431) we obtain 



^(^^\t)^^*{t + ta) + {l-^*{ta)) 



^^{t)-^^{t + ta) 



l-*?(ia) 




0.001 



ta = "= + 

t„ = 200 X 

t, = 50 M 

l, = 20 □ 

'a = 5 ■ 

'a = 2 O 

t = • 



• 8s u. O: 
• 8... » 

• a, u 
• a, ■: 

■■• 8. ■ 
• a^ 



1 



100 



(54) 



FIG. 1: The ta-old correlation function, $^''(t), for different 
values of ta- Curves represent a numerical integration of II40II 
witli 2 — 5/3 and T — 1.5, while dots correspond to the 
approximate formula II54II . 



V. CONCLUDING REMARKS 

The central result of this paper is the generalization 
of the GME discussed in Ref. [rj to one of arbitrary 
age. Another result, closely related to the GME of arbi- 
trary age, is the generalization of the Onsager principle 
to physical conditions of any age. The validity of this 
generalization of the Onsager principle is confirmed by 
the fact that Eq. 1(19(1 . the generalized Onsager princi- 
ple, yields Eq. 1(181) . and this is shown to be equivalent 
to the exact prescription of Godreche and Luck, which is 
independently re-derived in Section IIC. Another inter- 
esting result is given by Eq. ((54(1 . This is an analytical 
expression for ta-old correlation function, whose accuracy 
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has been established using Eq. (|40|l which is another ex- 
act expression for the ta-o\d correlation function, deter- 
mined by the Liouville-like approach of Section III. In the 
special case of non-Poisson processes where the function 
ip*{t) has the form of Eq. H32() . the approximate expres- 
sion Eq. (|54|l turns out to be very accurate. Of course 
there might be non-Poisson processes where Eq. (|54|) is 
not as accurate as in the case presented here. However, 
the iterative procedure discussed in Section IVC, allows 
us to determine higher-order corrections, should they be 
necessary for a more satisfactory treatment. 

It is important to understand why Eq. (|54|l is not ex- 
act, in spite of the fact that the approximation made for 
its derivation seems to fit the renewal nature of the pro- 
cess under study, where any jump resets memory to zero. 
This is a consequence of the infinite memory generated 
by non-Poisson dynamics, in spite of the fact the ran- 
dom events reset the memory to zero. The evaluation of 
the correlation function involves probabilistic arguments, 
and with them the infinite memory associated with the 
probabilistic treatment of non-Poisson processes. 

We illuminate the meaning of the correlation function 
of arbitrary age by means of the Liouville-like density ap- 
proach. The observation of the process of regression to 
equilibrium of the population difference corresponds to 
evaluating the anti-symmetric distribution, while leav- 
ing the symmetric part of the distribution free to evolve. 
If the distribution is at equilibrium, the symmetric part 
corresponds to the equilibrium distribution and the inte- 
gral of the left portion of the anti-symmetric part, with- 
out back injection, regresses to zero as the corresponding 
equilibrium correlation function. For any other condi- 
tion, the integral of the left portion of the anti-symmetric 
part regresses to zero with an analytical expression de- 
pending on the time at which observation begins. The 
regression continues as a function of that specific initial 
condition while the symmetric part keeps moving towards 
equilibrium independently of the population difference. 
This explains why the regression to equilibrium depends 
on the initial condition, of any age, with no further de- 
pendence on the bath dynamics that keeps drifting to- 
wards equilibrium. This also explains why an emission 
or absorption spectrum |23j | is not stationary and changes 
with time. The resonant radiation establishes a connec- 
tion between the anti-symmetric and the symmetric parts 
of the distribution, thereby updating observation to the 
changing bath conditions. 

It is worth ending this paper with some further remarks 
about these theoretical problems. We have built up a 
GME of arbitrary age, using an empirical approach. Is it 
possible to derive the same GME by using a Liouville-like 



approach? In principle, we should use the Liouville-like 
picture of Section III, to derive, via contraction on the 
bath variables, the same GME, of arbitrary age, as that 
of Section II D. However, it is evident that this effort, 
even if we were successful, would be of limited help, for 
practical purposes. Suppose, for instance, that we have 
to study the response of the system to an external, time 
dependent, perturbation. Would the GME of arbitrary, 
but fixed, age, useful for this purpose? It is evident that 
it would not. In fact, the external perturbation at times 
different from the age of the system would produce effects 
departing from the more realistic approach resting on 
perturbing trajectories. In the specific case of the absorp- 
tion spectrum of blinking quantum dots the authors 
adopted in fact this trajectory perspective to make a the- 
oretical prediction that is incompatible with the pertur- 
bation of a GME of fixed age. In literature, there already 
exists at least the discussion of one case [16] that seems 
to be a natural consequence of this property. Sokolov, 
Blumen and Klafter [l^ derived an exact density equa- 
tion to describe a sub-diffusion process. This equation 
corresponds to brand new initial conditions. This means 
a condition where ta = 0. Then, these authors perturbed 
this equation with a time dependent field, and they found 
that the theoretical result conflicts with the behavior of 
the CTRW under the influence of the same perturbation. 
This is so because the time dependent perturbation corre- 
sponds to additional observation, taking place at different 
time values, none of them coinciding with the observation 
time, but the perturbation at i = 0. This is true, what- 
ever the observation time is, either ta — 0, as in Ref. [l6j |. 
or cx) > > 0, a condition requiring the GME of this 
paper. Regardless of the observation time that we assign 
to the GME, it is impossible to make the GME predic- 
tion identical to the CTRW prediction, if we require the 
perturbation to remain external to the system. The con- 
cept of perturbation itself turns out to be inadequate to 
study non-Poisson processes, regardless of its intensity. 
Thus, the results of this paper, in addition to shedding 
light into those of Ref. [ig , imply a violation of the lin- 
ear response theory. The only possible way to make the 
density compatible with the trajectory picture is to make 
the external perturbation become a part of the system to 
study. This means that we have to build up a totally new, 
field-dependent, GME, along the hues of Ref. ^2^. This 
sets a limit on the applicability of the GME of arbitrary 
age found in this paper. However, this result seems to 
support the conclusion that the trajectory-density con- 
flict, revealed by Bologna, Grigolini and West !2l| might 
be a consequence of the aging properties emerging from 
non-Poisson renewal process. 
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